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SUMMARY  AND  CONCLUSIONS 
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1.0 


This  is  the  final  TSC  report  on  our  overview  of  the  TAGSEA  sea  clutter 
measurement  program.  It  includes  some  general  comments  on  the  techniques 
used  by  Raytheon  and  General  Dynamics  in  obtaining  the  data,  the  data  reduction, 
and  test  results.  Since  the  early  discussions  on  procedures  and  tests  were 
covered  in  monthly  meetings  and  reported  earlier,  this  report  emphasizes  the 
validity  of  the  results  and  their  use  in  relating  sea  clutter  model  for  Navy 


use. 


,/C.  ) ' 


t, 


TSC  has  no  rsason  to  question  the  validity  of  the  data  Itself.  While 
the  system  was  not  dssigned  for  quantitative  measurement  of  aQ , the  results 
should  be  accurate  to  within  a few  decibels  t£3  dB(2o)].  Both  Raytheon  and 
Ganaral  Dynamics  appear  to  have  performed  good] teats  in  the  time  allotted.  Two 

ara^s  deserve  special  attention.  1 V .. 

. -i  , '•  - 7 
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1.1  Conclusions  on  Reflectivity 

f'4 

While  the  measurement  of  sea  reflectivity,  (S' , was  not  a primary 

i ! y 

function  of  the  tests,  the  results  showed  the  greatest  deviation  from  earlier 

models.  For  grazing  angles  of  10,<^<45,l  the  TAGSEA  mean  values  of  crQ  were 

6-9  dB  higher  than  earlier  models,  especially  those  of  NRI  for  the  higher  sea 

states  (SS>4) . It  is  proposed  in  Section  3 that  TSC  believes  that  the  TAGSEA 

results  are  probably  more  representative  of  high  sea  state  reflectivity  than  the 

earlier  models  and  that  some  values  in  new  "models"  should  be  raised  by  4-6  dB.  v 

This  conclusion  is  not  solely  based  on  TAGSEA  results,  but  also  on  other  recent 

experiments  by  several  capable  experimenters. 


Aa  expected,  the  reflectivity  was  higher  on  vertical,  than  horizontal 
polarization.  Reflectivity  was  higher  in  the  upwind  direction  than  in  the  cross 
wind  direction  as  in  earlier  tests.  However,  the  downwind  reflectivity  data 
was  about  equal  to  the  upwind  TAGSEA  data.  Some  earlier  data  indicate  larger 
upwind/downwind  ratios  for  horizontal  polarization. 
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1.2 


Conclusions  on  Statistics 


The  TAGSEA  statistical  fluctuation  model  appears  to  be  a good 
engineering  design  model.  There,  the  total  fluctuation  ie  given  as  Che  product 
of  three  processes: 

0 ■ o «F*S 

O 0 

where  oq  is  the  average  cross-section 

S is  a slowly  fluctuating  process  representing 
long-term,  large-spatial  fluctuations  in  the 
■ cross-section 

and  F Is  a fast  fluctuation  process  which  has  a 

non-Rayleigh  distribution. 

The  above  model  is  consistent  with  a time-varying  Rayleigh  (TVR)  fluctuation 
model  which  has  been  proposed  by  several  authors  and  developed  in  Borne  detail 
by  Sodergren  [8]  and  Trunk  [22].  In  the  TVR  model,  the  primary  scattering 
mechanism  is  Bragg  scattering  from  (small  scale)  capillary  waves.  The  variable 
nature  of  the  slow  fluctuations  arises  from  local  tilt  produced  by  the  (large 
scale)  gravity  waves  and  by  variations  in  the  local  wind  field.  A usual  appli- 
cation of  the  TVR  distribution  models  the  fast  fluctuation  as  a Rayleigh  process. 
Such  application,  however,  ignores  the  difficulty  in  obtaining  Independent  samples 
from  the  same  range  cell  under  "frozen" conditions.  That  is,  under  conditions 
such  that  the  local  tilt  or  wind  field  does  not  change  during  the  sampling 
interval.  The  TAGSEA  model  provides  a hedge  against  this  mistake  by  lumping 
some  of  the  non-Rayleigh  characteristics  into  the  fast  fluctuations. 

As  noted  in  [1],  the  model  is  limited  to  a fairly  specific  set  of 

conditions  — high-resolution  radar,  high  sea  state,  and  moderate  grazing  angle. 

Comparison  with  other  data  indicates  that  at  low  grazing  angle  the  o fluctuations 

-5  ° 

have  much  higher  tails  (2-8  dB  at  10  probability  of  a threshold  crossing)  and 
show  fairly  strong  dependence  on  sea  state.  The  TAGSEA  model  should  only  be 
applied  to  conditions  similar  to  the  data  set  from  which  it  was  derived.  Compar- 
ison with  the  NRL  20  nsec  radar  indicate  that  a higher  resolution  radar  may  produce 
higher  tails  than  the  TAGSEA  model.  Thus,  a conservative  design  practice  would 
be  to  design  to  a somewhat  more  severe  model  (e.g.,  1-2  dB  higher  tails  at 
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Theoretical  calculations,  Sodorgren's  model,  and  the  TAGSEA  model 
differ  on  the  rate  of  rapid  fluctuations.  For  coherent  frame  rates  that 
are  of  the  same  order  aa  the  fast  decorrelation  time,  both  the  amplitude 
distribution  of  clutter  signals  (which  affect  detection  thresholds)  and 
second  order  statistics  (which  affect  MTl/Doppler  processing)  would  depend 
on  the  integration  and  decorrelation  processes. 


I 

1 2.0  TECHNOLOGY  SERVICE  CORPORATION  TASKS 

| TSC ’a  role  in  Che  TAGSEA  program  was  divided  into  two  areas.  The 

• fir at  role  occurred  during  the  measurement  and  algorithm  development  phases. 

, Xt  consisted  of  reviewing  the  objectives  and  plans  in  both  of  these  areas.  The 

! mode  of  operation  was  primarily  by  direct  technical  discussions  on  a monthly 

basis.  These  are  summarized  in  minutes  of  those  meetings. 

The  second  role  occurred  both  during  and  after  the  measurements 
I and  consisted  of  reviewing  the  data  as  available  from  Raytheon  and  supplying  them 

with  existing  models  for  them  to  check  their  initio.’  results.  The  review  of  the 
Raytheon  final  report  is  given  in  Section  3.0  of  this  report. 

The  TSC  tasks  are  shown  below: 

\.  "TSC  will-  review  the  Raytheon  Test  Objectives  and  Test  Plan  to  Insure  

that  the  data  obtained  will  contain  adequate  information  to  determine 
the  desired  sea  clutter  statistics.  This  will  Include  TSC  evaluations 
of  the  PRF  and  the  number  of  independent  samples,  dynamic  range  optimi- 
zation, allowance  for  aircraft  motion,  polarization  properties,  sufficiency 
of  the  number  of  samples  and  sea  conditions  within  the  constraints  of  the 
overall  program. 

"TSC  will  review  the  instrumentation  and  especially  the  calibration 
techniques.  TSC  will  help  insure  that  the  critical  problem  of  obtaining 
sufficient  sea  state  descriptors  are  in  the  planning  at  the  time  of  the 

tests. 

"TSC  will  review  the  data  reduction  techniques  in  the  critical  areas 
of  obtaining  distribution  functions  for  both  small  areas  resulting  from 
beam  sharpening  in  azimuth,  and  the  larger  regions  defined  by  the 
azimuth  beamwidth. 

B.  "TSC  will  Interpret  the  results  in  the  light  of  Navy  requirements. 

Verification  of  the  results  will  be  performed  both  by  reviewing  the 
procedures,  data  reduction  and  analysis  described  in  Tasks  1 and  2 and 
by  comparing  certain  results  with  earlier  reflectivity  models  developed 
by  NRL,  TSC,  APL  and  others.  In  addition,  the  measured  distributions 
will  be  checked  against  models  developed  by  NRL  (Trunk) , APL  (Sodergren) 
and  in  unpublished  TSC  reports." 
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A Final  Report  will  be  issued  that  includes  verification  of  the 
Raytheon  results  and  comparison  with  other  models. 

ISC  personnel  attended  meetings  at  Raytheon  on 

12  January  1976 
25  February  1976 
12  April  1976 
11  May  1976 
22  June  1976 
30  July  1976 

plus  several  meetings  at  AFL/JHU.  TSC  memos  related  to  these  meetings  are 
available  [ 25,  29,  30,  31,  32,  33]. 

2.1  Validation  of  Teat  and  Data  Reduction  Procedures 

In  general,  the  constant  interaction  between  APL/JHU,  General 
Dynamics,  TSC  and  Raytheon,  Bedford  achieved  a set  of  teat  procedures  and  data 
reduction  techniques  that  were  acceptable  to  all  parties.  It  must  be  pointed 
out  that  the  distributions  presented  are  describing  a non-stationary  process, 
and  there  is  no  established  "right"  or  "wrong"  in  the  selection  of  observation 
times,  number  of  samples,  and  spatial  averaging. 

TSC  would  have  preferred  slightly  more  emphases  on  distributions 
of  a single  spatial  cell,  but  this  does  not  seem  to  have  affected  the 
conclusions. 

TSC  would  like  to  have  obtained  more  data  on  horizontal  polarization 
and  at  lower  grazing  angles,  but  the  short  time  scale  and  equipment  limitations 
did  not  permit  this. 


3.0 


VALIDATION  AND ' INTERPRETATION  OF  RESULTS 


This  section  relates  the  reported  Raytheon  data  to  both  prior 
experiments  and  theory.  This  section  is  split  arbitrarily  between  mean 
reflectivity,  oq  (or  o°),  and  sea  clutter  statistics.  However,  a revised 
theoretical  model  for  reflectivity  is  Included  in  the  statistics  section. 


f • 

1 

1 

I 

I 

I 

I 


3* 1 Reflectivity  of  Sea  Clutter 


While  sea  clutter  reflectivity  measurements  have  been  performed  for 
over  30  years,  controversy  remains  in  both  the  reflectivity  and  statistical 
models  in  some  specific  areas.  The  general  rules  and  their  limitations  for 
grazing  angles  of  l*-40*  can  be  summarized  as  follows  for  a fully  developed 

sea: 


....a.  The  monostatic  b-ckscatter  coefficient,  oq,  monotonically  j 

increases  with  grazing  angle,  i K from  below  1°  to  above  .j 

40*  for  all  sea  states,  with  a higher  slope  from  1°-10* 
and  a lesser  slope  from  10*-40*  (depending  on  frequency  and  sea  state). 


b.  The  backacatter  coefficient  monotonically  increases  with 

sea  state  with  a higher  slope  at  lower  grazing  angles, 

lower  carrier  frequencies,  and  lower  sea  states.  The 

rate  of  change  of  oQ  decreases  at  the  higher  sea  states. 

There  have  been  claims  that  at  the  higher  frequencies 

(f  >5GHz)  there  is  a "saturation"  at  about  Sea  State  4. 
o 

This  is  the  subject  of  active  controversy  between  NHL 
(26,  34  ] which  believes  in  saturation  and  NASA  related 
experlmentors  and  theoreticians  [ 5,  12  ] who  are  estimating 
wind  speed  (and  sea  state)  from  reflectivity  for  satellite 
radar  programs.  The  NASA  experiments  and  related  theory  by 
Pierson  show  no  saturation  up  to  about  Sea  State  7. 

- c.  The  reflectivity  increases  monotonically  with  carrier 

frequency,  fQ,  from  0.5  GHz  to  at  least  5.0  GHz.  It  also 

appears  that  o decreases  with  frequency  from  f <35  GHz 
o o 

to  over  fQ“100  GHz.  At  this  point  there  is  again  controversy 
with  NRL  measurements  at  high  sea  states  tending  to  decrease 
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at  10  GHz  from  measurements  at  5GHz  while  other  experi- 
ments show  Increasing  with  frequency  to  above  10  GHz 
and  perhaps  to  15-18  GHz  or  even  35  GHz. 

d.  The  mean  reflectivity  on  vertical  polarization  exceeds  that 
of  horizontal  polarization  for  1®<^<40°,  and  virtually  all 
sea  states  for  carrier  frequencies  of  0.5  to  about  15  MHz. 

The  amplitude  distribution  is  also  closer  to  a Rayleigh 
distribution  for  vertical  polarization  than  is  the  more 
"spiky"  horizontal  polarization.  ThuB,  a caution  must  be 
attached  to  using  the  lower  values  of  cq  on  horizontal 
polarization  when  automatic  detection  circuitry  is 

contemplated.  False  alarm  problems  may  develop  when 
attempting  to  detect  small  aircraft  or  ships  in  a spiky 

sea  clutter  background.  Fig.  3.1  illustrates  the  effeet  using 

the  Raytheon  "TAGSEA  Extreme  Model".  The  ordinate  is  the 

probability  that  aQ  exceeds  the  values  on  the  abscissa. 

The  Important  point  is  that  while  in  thiB  example,  the 

median  value  of  a is  6 dB  lower  for  c*  than  for  a*  , the 
O Hli  W 

receiver  threshold  level  required  to  limit  the  clutter 
false  alarm  probability  to  10“5  differs  by  only  2 dB. 

Further  information  on  thaise  statistics  is  given  In  Section  3.2 

e.  The  reflectivity  is  higher  when  looking  upwind  rather  than 
croaswlnd . These  values  are  typically  2-6  dB.  Downwind 
values  are  generally  between  the  two,  but  closer  to  the 
upwind  value.  NASA  and  NRL  reports  have  considerable 
measured  data  with  some  in  a recent  Boeing  Report  [9]. 

Note  that  all  of  the  above  refer  to  the  grazing  angle  of  the  EM 
radiation  and  not  the  depression  angle  of  the  radar  aperture.  Note  also  that 
the  data  is  primarily  for  pulse  durations,  t,  greater  than  0.25  microsecond 
or  125  ft  in  radial  extent.  Effects  of  beamwidth  or  azimuthal  extent  are 
suspected  but  not  found  to  date. 

Without  TAGSEA  data,  there  waB  controversy  about  the  mean  values  and 
the  statistics  in  certain  situations.  The  TAGSEA  experiments  provided 
additional  data.  These  are  presented  in  the  next  section. 
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3.1.1 


TAGSEA  Results  and  Data  Comparison 

The  TAGSEA  data  on  reflectivity  can  be  summarized  quite  easily: 


frequency  ■ X-Band 


grazing  angles 
polarization  ■ 

resolution  ■ 

aea  state  • 


4. 7®<i|i<52.8* 

both  linear  polarizations  (more 
vertical  data  than  horizontal) 

118  ft  effective  in  range  and  varying 
in  cross  range,  but  typically  100' 

1 to  5 + with  reasonable  receiver 
•ignal-to-noise  ratios  at  the  higher 

aea  states 


The  results  for  <J  were  self  consistent  from  run  to  run,  and  for 
0 

measurements  in  the  Atlantic  and  Pacific  Oceans.  The  values  of  a were  consist 

o 

tenjly  6-10  dB  higher  than  the  NRL  1965-1969  data  which  forms  most  of  earlier 
data  base. 


The  summarized  Raytheon  TAGSEA  results  from  Fig.  4-8  of  BR-9254-1 
Volume  I (1]  are  shown  on  Fig.  3.1-1  along  with  the  earlier  TSC  projections  for 
the  same  conditions.  Note  the  consistently  higher  reflectivity  values  from  the 
TAGSEA  tests. 

It  must  be  pointed  out  that  the  earlier  TSC  models  are  to  a great  extent 
the  NRL  data  taken  from  the  4FR  aircraft  in  the  1960's.  The  number  of  data 
points  from  this  aircraft  was  so  great  that  it  far  outweighs  other  data  points  at 
X-Band  carrier  frequencies. 

Some  of  the  NRL  data  points  for  '('■10®  \pm309  are  illustrated  on  Figures 
3.1-2  and  3.1-3.  While  the  concept  of  a "saturation"  in  cjo  as  sea  state  increased 
to  about  SS4,  bothered  a lot  of  theoreticians,  the  reliability  of  the  equipment  and 
general  care  in  calibration  led  people  to  chwose  the  experimentally  derived 
models. 


However,  in  the  early  1970's,  NASA  was  engaged  in  a program  to  measure 
wind  speed  by  sea  reflectivity  measurements  using  a "scatterometer",  a radar 
technique  to  estimate  These  were  part  of  the  SKYIAB  and  SEASAT  programs. 

If  the  saturation  theory  was  correct,  their  device  could  not  predict  high  winds 
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(or  high  sea  stai.es).  NASA,  with  the  University  of  Kansas,  flew  many  flights 
of  a Ku“Band  radar  (13.9  GHz)  to  measure  ao  vs  wave  height  and  wind  speed. 

They  also  engaged  Dr.  Willard  Pierson  of  NYU  to  relate  radar  echoes  to  ocean- 
ographic parameters.  The  NASA  data  and  related  theories  did  not  exhibit  the 
saturation.  Thus,  two  extensive  sets  of  data  existed  with  the  only  difference 
being  the  carrier  frequency  difference  between  9.5  and  14  GHz. 

In  order  to  validate  the  TAGSEA  results,  the  NASA  data,  along  with 
the  NRL  and  TAGSEA  results  and  other  recent  measurements,  are  graphed  on 
Figures  3.1-4  and  3.1-5  for  Horizontal  and  Vertical  polarization.  The  primary 
difference  with  other  earlier  graphs  is  that  high  sea  state  data  points  at 
several  carrier  frequencies  are  shown  on  the  same  graph/  Also  Included  are 
APL  data  pointB  [8,  35]  that  had  not  heretofore  been  available  to  TSC,  and 
recent  S-Band  measurements  by  Boeing  Aircraft  Corporation  [9].  Note  that  the 
X-  and  Ku-Band  data  show  little  difference  except  that  now  the  NRL  oQ  data 
seems  low.  The  10  dB  lower  values  for  the  Boeing  S-Band  data  can  probably  be 
explained  by  the  factor  of  3 difference  in  carrier  frequency. 

Thus,  the  saturation  effect  of  the  NRL  data  appears  to  be  violated 

by  several  experiments  including  TAGSEA.  Since  'the  NASA  data  set,  which  is 

the  second  most  extensive,  was  acquired  at  a somewhat  higher  carrier  frequency 

and  only  overlaps  the  NRL  data  at  the  higher  grazing  angles,  it  alone  would  not 

be  sufficient  to  alter  the  models.  However,  the  additional  APL  data  at  both 

X-  and  Ku~Bands  and  the  Sittrop  data  [28]  seem  to  confirm  the  higher  values  for 

o.  at  high  sea  states  and  that  there  is  little  difference  in  c between  X-  and 
o o 

K^-Band . APL  72  refers  to  references  [23]  and  [35]. 

Some  care  must  be  used  in  interpretation  since  cq  is  generally  given 
as  a mean  value  while  NRL  generally  publishes  a median  value.  The  Boeing  data 
1b  also  a median  value.  The  mean- to-med lan  ratio  is  typically  +1.6  to  2.0  dB, 
for  the  angles  of  interest,  and  perhaps  higher  at  the  lower  grazing  angles. 

This  has  the  effect  of  raising  the  NRL  and  Boeing  data  and  closing  the  gap. 

A model  is  ulscussed  in  the  next  sections. 

*The  winds  were  all  in  excess  of  22  knots  with  the  exception  of  the  APL  Tests. 
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3.2 


Evaluation  of  the  TAGSEA  Model  of  Saa  Clutter  Statistics 


In  this  section,  the  validity  of  the  statistical  fluctuation  model 
proposed  in  reference  [ll  is  examined.  Although  absolute  validity  is  extremely 
difficult  to  establish,  it  is  possible  to  examine  the  correspondence  between 
the  model  and  existing  information  about  the  reflectivity  of  the  sea  surface. 
Generally,  a model  should  be  consistent  with  accepted  theory  of  backscatter 
from  a rough  surface,  and  should  agree  with  experimental  measurements.  If 
the  model  satisfies  these  criteria,  then  until  an  improved  theory  is  developed, 
or  new  experimental  observations  are  made,  it  is  usually  accepted  as  a valid 
characterization.  In  the  following  sections,  these  aspects  of  the  TAGSEA 
fluctuation  model  are  examined. 

Generally  it  is  found  that  the  TAGSEA  model  agrees  well  qualitatively 
with  1 time  and  space  varying  Rayleigh  distribution.  Such  a model  can  be 
derived  from  a composite  sea  surface  model  for  which  gravity  waves  tilt  the 
local  surface  seen  by  the  radar  and  for  which  the  primary  scattering  mechanism 
la  Bragg  scattering  from  capillary  waves.  In  addition  to  qualitative  agreement, 
significant  quantitative  agreement  exists  for  vertical  polarized  cases  when 
the  Pierson  and  Stacy  Cs]  sea  surface  model  is  used.  Since  that  model  was 
derived  largely  from  direct  water-wave  measurements  rather  than  radar  measure- 
ments, the  agreement  tends  to  vslidets  both  models. 

A comparison  with  experimental  data  lndlcatai  that  most  previously 
reported  data  has  been  taken  at  low  grazing  angles  and  ia  not  appropriate  for 
comparison  with  the  TAGSEA  model.  Four  references,  [6]-[9]  were  available 
to  provide  high-resolution,  moderate-grazing-angle,  and  hlgh-sea-state  data 
for  comparison  with  the  TAGSEA  model.  Of  these  data,  measurements  from  the 
NRL  20  nsec  radar  gave  distribution  functions  for  tha  cross  section  fluctuation 
which  had  significantly  higher  tails  than  the  model.  The  Boeing  data  was  too 
preliminary  for  detailed  comparison  but  appeared  to  give  results  similar  to 
the  TAGSEA  model.  Sodergren's  results  provided  the  closest  agreement  with 
the  model. 


Generally  the  model  proposed  in  [l]  appears  to  be  good  for  engineering 
design  and  evaluation.  As  was  noted  in  [l],  it  was  derived  for  a specific  set 
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of  condition*  — high-raaolutlon-radar,  moderate-grazing-angle,  and  high-aaa- 
atate.  Comparison  with  data  taken  at  low  grazing  angle  [8],  [lO],  [ll] 
indicates  that  significant  modification  of  the  model  is  required  at  low 
grazing  angle,  In  this  caso,  the  model  for  fluctuation  statistics  will 
be  sea  state  dependent,  and  will  generally  have  significantly  higher  tails 
than  the  TAGSGA  model  (e.g.,  as  much  as  2-8  dB  at  the  10  ^ point,  depending 
on  sea  state). 

No  consistency  has  been  found  for  determining  the  rapidity  of  the 
fast  fluctuations,  The  TAGSGA  model,  Sodergren's  modal  and  the  theoretical 
discussion  of  Section  3.2.1  all  give  different  relations  for  the  bandwidth 
of  the  fast  fluctuations.  For  low  data  rata  systems  (frame  rate  <100  Hz) 
this  is  not  Important  since  all  the  above  models  degenerate  to  independent 
fluctuations.  For  very  high  data  rate  systems,  the  problem  will  become 
significant  and  should  be  investigated. 

3.2.1  Comparison  Between  the  TAGSGA  Modal  and  Rough  Surface  Scattering 

Theory 

Theoretical  models  for  scattering  from  rough  surfaces  have  been 
studied  by  many  investigators  both  from  the  atandpoint  of  electromagnetic 
scattering  from  above  the  sea  surface  and  acoustic  scattering  from  below. 

The  references  on  the  subject  are  extensive,  and  those  listed  [l2]-[l6] 
represent  only  a small  sample.  In  the  reference  list,  Swift  and  Jones  [12] 
give  a nice,  simple  discussion  of  electromagnetic  scattering  from  the  see 
surface.  Peake  and  Oliver  [l3],  Barrlck  [14],  or  Barrick  and  Peake  [15] 
have  more  comprehensive  discussions.  Horton  end  Muir  [16]  give  e good 
summary  of  the  theory  of  acoustic  scattering  from  the  aaa. 

The  theoratlcal  model  that  appears  to  be  gaining  the  widest 
acceptance,  especially  in  remote  sensing  and  satellite-radar-scatterometry , 
is  the  composite  surface  model  (see  Peake  end  Oliver  [l3],  pg.  74-77).  In 
this  model,  the  ocean  surface  is  characterized  as  the  sum  of  two  Independent 
processes  - a large  scale  roughness,  which  has  a small  scale  roughneas 
superimposed,  Wright  [24]  has  applied  this  model  to  a series  of  radar 
measurements.  For  scattering  at  X and  band,  the  large  acale  roughness 
is  approximately  associated  with  the  gravity  waves  and  swell,  while  the 
small  acale  roughnaaa  ia  associated  with  the  capillary  wavos. 
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Generally  "large  scale  roughness"  Is  defined  to  mean  tlint  tha 
surface  height  variations  are  large  relative  to  the  rndnr  wavelength. 
Scattering  from  this  scale  can  be  computed  when  the  surface  is  gently 
undulating  so  that  the  curvature  is  always  much  greater  then  the  radar 
wavelength,  and  the  slopes  of  the  surface  are  small.  When  these  conditions 
hold,  the  scattering  properties  can  be  obtained  by  computing  the  number  and 
curvature  of  the  specular  points  as  is  illustrated  in  Fig.  3.2-1  (Peake  and 
Oliver  [l3]).  Barrick  [14]  has  analyzed  this  case  in  some  detail,  and  has 
derived  mean  cross  section  as  well  as  cross  section  statistics. 

For  a rough  surface,  the  mean  cross  section  is  given  by 


where 


ao  _ .[Rl^pc4  9 exp  [..((tan  e)/S)2]  (1) 

S* 

6 is  the  angle  of  incidence,  i.e.,  the 
angle  between  the  look  direction  and 
the  normal  to  the  surface 


and 


f\ 

|r|*  is  the  power  reflection  coefficient 
of  a plane  wave  normally  incident  on 
a flat  dielectric  surface  (at  X band 
|R|2  is  approximately  -2.1  dB  for  sea 
water) 

2 

S is  the  rms  slope  of  the  surface. 


The  normalized  cross  section  versus  Incident  angle  is  given  in  Fig,  3.2-2  for 
a Gaussian  surface  as  a function  of  the  rms  slope. 


For  a fully  developed  sea  surface,  Swift  and  Jones  [12]  give  the 
rms  slope  as  a function  of  wind  velocity  in  meters  per  second  as: 


S 


2 


(3.16  x 10  "*)  V (up/down  wind) 
.003  + .00192  V (cross  wind) 


(la) 
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Figure  3.2-1.  Scattering  from  an  Undulating  Rough 
Surface;  h » X,  slopes  « 1,  scattered  beamwidth 
m surface  slope,  specular,  (from  [l3]) 
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Figure  3. 
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2-2.  Average  Radar  Croat  Section  per  Unit  Surface  Area 
0°  for  Very  Rough  Gouaaian  Surface,  (from  [14]) 
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Note  that  for  sea  state  5,  a typical  wind  velocity  la  about  11  m/a  (22  knots); 
this  gives  an  rms  slope  of  about  10°.  From  Figure  3.2-2  one  would  expect  the 
rough  surface  scattering  cross  section  to  be  small  at  incidence  angles  greater 
than  30  degrees.  All  the  TAGSEA  results  correspond  to  higher  angles  of  incidence; 
thus  the  specular  scattering  model  which  applies  near  vertical  incidence  would 
not  be  expected  to  apply  to  the  TAGSEA  results. 

The  small  scale  roughness  is  defined  to  be  those  height  variations 
which  are  much  smaller  than  the  rf  wavelength.  The  scattering  mechanism  from 
this  component  of  the  surface  is  completely  different  than  the  rough  surface 
scattering  discussed  previously.  The  case  of  a slightly  rough  surface,  can  be 
computed  by  a perturbation  technique  described  by  S.O.  Rice  [17],  the  first 
order  coefficients  have  been  developed  in  detail  by  Barrick  [id], 

.In  general,  the  total  field  produced  by  e slightly  rough  surface 
is  the  sum  of  a coherently  reflected  and  a diffusely  scattered  component  as 
sketched  in  Fig.  3.2-3.  The  coherent  component  is  identical  to  the  field 
reflected  from  a perfectly  flat  eurfeca,  but  the  power  has  been  reduced  by 
the  factor  axp(-2(khcos6i)'i)  where  k is  the  rf  wave  number,  and  h la  the  rm* 
surface  height,  The  first  order  exprassion  for  the  diffuse  component  can  be 
shown  to  be  equivalent  to  the  Fourier  transform  of  the  surface  height. 

A useful  concept  for  understanding  the  characteristics  of  the 
diffuse  component  is  to  consider  each  point  on  the  slightly  rough  surface  as 
a nearly  isotropic  scattering  point.  The  field  which  is  reradiated  is 
proportional  to  the  incident  field,  but  the  phase  is  determined  by  the  height 
of  the  surface.  Such  a surface  acta  as  a "phased  array"  with  small  random 
phase  errors  which  are  proportional  to  the  height  of  the  surface.  Since  the 
array  is  "fed"  by  the  incident  plane  wave,  it  la  phased  so  that  the  maximum 
response  is  in  the  specular  direction.  The  diffuse  component  corresponds  to 
sidelobes  produced  by  the  random  phase  errors.  For  small  enough  phase  errors 
(i.e.,  a slightly  rough  surface),  the  sidelobe  response  is  proportional  to  the 
Fourier  transform  of  those  errors,  Note  that  the  diffuse  component  is  not  the 
result  of  many  Independent:  point  scatteres,  Bince  the  phase  and  amplitude  of 
each  point  on  the  surface  is  determined  by  the  Incident  and  look  directions, 


3-16 


and  the  height  of  the  surface.  Although  the  Burface  is  random,  it  has  a finite 
correlation  distance,  so  that  the  scattering  points  are  not  independent. 

Significant  scattering  power  occurs  in  those  directions  which  have 
significant  Fourier  components  in  the  surface.  This  is  usually  known  ns  Bragg 
scattering.  In  general , if  the  surface  is  composed  of  components,  the  scatter- 
ing is  proportional  to  those  components  which  resonate  with  the  vector  difference 
between  the  propagation  directions  of  the  incident  and  scattered  plane  wave 
numbers . 


For  backscatter  Barrick  [l4]  gives  the  average  normalized  ctobs 
section  for  a slightly  rough  surface  as: 

0°  * 4nk^cos^6,^  |»  | (»2ksin0^,O)  (2) 

where 

2tt 

k « T“  , is  the  RF  wave  number, 

0^  is  the  incidence  angle  (measured  from  the  vertical), 

a is  a scattering  coefficient  which  is  a function  of 
the  surface  materials,  polarization  states  and 
incidence  angle. 

and 

Wft.m)  is  the  two-dimensional  Fourier  transform  of  the  two- 
dimensional  autocorrelation  function  of  the  surface, 

The  scattering  coefficient,  O',  for  a nonmagnetic  surface  is: 
e-1 

O'.  . ■ — - -g. 

(cos  0^  + Jc  - sin^  0^) 

and 

* « (<!-!)[ (e -1) sin2  Ei  + e] 

vv  | r-  2 

(c  cos  0^  +s]  e - sin  0^) 

whore  c is  the  relative  permittivity  of  sea  water.  In  the  above  expressions, 
"hh"  and  "vv"  refer  to  horizontal  and  vertical  polarization  respectively. 
Table  I gives  e ns  a function  of  frequency  for  a nominal  water  temperature 
of  10°  centigrade. 
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Table  1.  Complex  Relative  Permittivity  of 
Sea  Water  at  10°C 

Frequency 

GHz 

Real  Part 

Imaginary 

.5 

83.6 

-144.0 

1.25 

82.9 

- 58.3 

3.0 

79.5 

- 28.2 

5.0 

73.1 

- 25.0 

9.7 

54.6 

- 36.5 

15.0 

37.9 

- 50.6 

35.0 

14.5 

- 71.3 

The  general  form  of  the  sea  cross  section  due  to  Bragg  scattering 
has  been  known  for  some  time.  The  greatest  unknown  in  equation  (2)  is  the 
magnitude  of  the  sea  surface  spatial  spectrum.  Pierson  and  Stacy  [5]  give 
a detailed  exposition  of  the  spatial  spectrum  of  the  sea  surface.  They 
identify  five  regions  of  surface  waves  covering  6 orders  of  magnitude  in 
wave  number  from  wavelengths  of  700  meters  to  .7  millimeters.  Using  avail- 
able data  from  direct  measurements  of  the  sea  surface,  photographs  of  waves 
generated  in  wind-watar  tunnels,  plus  theoretical  models,  Pierson  and  Stacy 
generate  a consistent  quantitative  model  of  the  sea  spectrum  for  each  of  the 
five  regions . The  resulting  model  is  dependent  only  on  the  wind  velocity 
over  the  sea  surface.**  Spectrums  of  the  one  dimensional  spectrum  from  the 
model  are  given  in  Figure  3.2-4  and  the  various  regions  are  marked.*  On 
that  figure,  U.Q  c represents  the  wind  speed  measure  at  19.5  m altitude,  and 

— • J 

Hj^  is  the  significant  waveheight.  Note  from  that  figure,  that  at  X and 
band,  the  Eragg  scattering  will  be  influenced  nriraarily  by  the  capillary  waves. 


* Since  we  are  concerned  primarily  with  the  capillary  waves,  only  his 
region  is  discussed.  The  intarested  reader  is  referred  to  [5]  for 
definitions  and  discussions  of  the  4 other  regions. 

**  For  a fully  developed  sea. 
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In  the  capillary  region,  Pierson  and  Stacy's  model  of  the  two 
dimensional  wave  spectrum  is  given  as* 

.*13  H 2 

W4  <<.,«) r-4'2  F(-f.,m)  (3) 

(<t  +m*') 


where 


and 


K and  m are  the  radar  wave  numbers  in  the 

upwind  and  cross  wind  direction. 

H is  the  significant  height  of  the  capillary 
c 

waves 

FOt,m)  gives  the  anisotropy  of  the  capillary  wave 
spectrum. 


In  Pierson  and  Stacy's  model,  the  capillary  wave  height  is  given  as 
a rather  complicated  function.  Figure  3.2-5  gives  versus  wind  speed. 

Although  it  is  relatively  safe  to  extrapolate,  the  heights  of  the 
capillary  waves  from  wind-water  tunnels  measurements  to  the  open  ocean,  it 
lv  not  possible  to  do  so  for  the  directional  characteristics  of  the  spectrum. 

The  problem  is  that  in  a water-wind  tunnel,  the  wind  direction  is  very  steady, 
and  the  resulting  capillary  waves  tend  to  be  generated  at  an  angle  to  the 
wind  direction  as  sketched  in  Figure  3.2-6.  This  results  in  a distinct  bimodal, 
directionality  to  the  wave  spectrum.  In  the  open  ocean,  the  wind  direction  is 
constantly  changing  and  tends  to  wash  out  the  bimodal  nature  of  the  spectrum. 

For  this  reason,  Pierson  and  Stacy  were  unable  to  determine  a unique  expression 
for  F(-l  ,«0  . 


To  obtain  a model  for  calculation,  this  author  decided  to  use  the 
direction  coefficients  inferred  by  Pierson  and  Stacy  from  some  NASA  upwind, 


**  In  aquation  (3),  the  subscript  4 is  for  region  4 which  is  the  capillary 

region.  The  reader  should  know  that  Burrlck  and  Pierson  differ  by  a factor 
of  2 in  the  definition  of  W.  The  relation  is 


W„  . , - 2Wnj  (-t.m) 

Barrick  Pierson  ’ 

the  expression  given  above  is  compatible  with  Barrick' s definition  so  that 
it  may  be  used  in  equation  (2). 
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croeawind,  radar  measurements  of  relative  c°  at  Ku  band  (Bradley  [l9]). 
Since  these  measurements  were  made  over  the  open  ocean,  at  a high  enough 
frequency  to  insure  capillary  wave  scattering,  it  was  felt  that  they  would 
be  most  suitable  for  comparison  with  the  TAGSEA  model.  From  these  measure* 
ments,  the  directivity  is  given  by 

F(f-,m)  ■ F(0)  - ~ [l  + a^  coa20]  (A) 


where 

0 

is  the  look  direction  relative 

to  the  wind  direction, 

and 

*1 

is  a coefficient  which  depends 
only  on  wind  speed, 

• « 4 

That  it, 

#i  " 

1.28  U*  23  -1 
1.28  IT23  +1 

where 

U 

is  the  wind  speed  at  19.5  meters 
height  in  knots. 

Substituting  (3)  and  (4),  into  the  expression  for  backscatter 
(eq.  (2))  yields: 

0°  - 0.167  |ff|2(cot4  61)(Hc2F(0))  (6) 

Note  in  the  above  expression  that  the  only  term  sensitive  to  RF  frequency  is 
|a|2.  Table  II  gives  |a|2  for  horizontal  and  vertical  polarization  as 
a function  of  grazing  angle  and  frequency.  That  table  show  negligible 
dependence  on  frequency;  thus  Bragg  scattering  from  capillary  waves  is 
dependent  only  on  polarization,  grazing  angle  and  wind  velocity. 

Equation  (6)  was  evaluated  for  conditions  comparable  to  the  TAGSEA 
measurements  and  the  results  are  given  in  Figs.  3.2*7  to  3.2*11.  The  comparison 
is  good  for  vertical  polarization  cases.  In  those  cases,  the  difference  between 
the  TAGSEA  measurements  and  the  theoretical  model  is  about  1*2  dB.  Ihe 
comparison  with  the  horizontal  polarization  measurements  is  not  as  good. 

The  primary  reason  for  this  is  that  the  model  has  not  yet  fully  included 
the  composite  nature  of  the  sea  surface.  In  particular  all  the  previous 
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Figure  3.2-8.  Comparison  of  TAGSEA 


Figure  3 . 2-11.  Coop arts an  of  TAGSEA  Data  with  Eq . s (6)  and  (IJ) 


derivations  wore  made  for  a slightly  rough;  flat,  horizontal  surface.  In  fact, 
the  capillary  waves  are  superimposed  on  a larger  gravity  wave  structure,  as 
sketched  in  Figure  3.2-12 

At  sea  state  5,  the  gravity  wavelength  is  on  the  order  of  60-80 
meters;  thus  for  the  short  TAGSEA  resolution  cell  (“  30  m) , the  portion  of 
the  gravity  wave  seen  is  fairly  flat,  but  it  is  not  horizontal.  Instead 
as  each  wava  passes  through  the  range  cell,  the  local  angle  of  incidence  varies 
according  to  the  slope  of  the  wave.  To  obtain  the  average  cross  sections,  the 
cross  section  computed  from  eq.  (6)  must  be  averaged  over  the  variations  in 
the  angle  of  incidence  produced  by  the  gravity  waves.  If  P(S)  is  the  probability 
density  distribution  of  the  wave  slopes,  the  average  crooB  section  is  given  by: 

w 

00cet)  - / O (81  - tan”1 (S) ) PCS)  dS  (7) 

0 

The  above  computation  will  have  little  effect  on  the  mean  cross  section 
for  vertical  polarization  since  C0^  varies  slowly  with  incidence  angle.  However, 
it  will  have  a significant  effect  on  the  horizontal  polarization  results. 

Eq.  (7)  was  evaluated  for  horizontal  polarization  assuming  a Gaussian 
slope  distribution  with  an  rms  slope  gain  be  eq.  (la).  The  results  are  plotted 
as  dotted  lines  in  Figures  3.2-10  and  3.2-11.  Note  that  the  agreement  is 
better,  but  still  not  as  good  as  the  vertical  polarization  results.  The 

reasons  for  this  are  not  known.  One  possibility  1b  the  horizontal  tilting  of 

the  gravity  waves.  That  is,  in  addition  to  tilting  "fore  and  aft",  the  gravity 
waves  also  tilt  "sido-to-eida" . When  the  surface  is  tilted  side-to-side , the 
incident  electric  field  is  no  longer  exactly  horizontally  polarized  with  respect 
to  the  local  surface  tilt.  Thus,  locally,  the  wave  has  slight  vertical 
polarization.  Since  vertically  polarized  backscatter  at  low  grazing  angle  is 
10-15  dB  higher  than  horizontal  polarization,  a slight  "local  vertical" 
component  will  significantly  change  the  <7°  characteristic.  Unfortunately  the 
computations  required  to  do  the  two  dimansionul  averaging  were  too  tedious 
to  complete  for  this  report.  It  should  be  noted  that  in  his  model,  Wright  C'24] 

did  include  tilting  in  both  directions.  He  also  was  able  to  obtain  good 

agreement  for  vertical  polarization  but  stated  that  tho  horizontal  polarization 
results  were  not  ns  good.  There  are  several  differences  in  the  details  between 
the  model  used  here,  and  Wright's  model,  but  tho  primary  difference  is  in  the 
assumed  water-wave  spectrum. 
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See  Surface 


Figure  3 


local  Tilt 


2~12.  Local  Tilting  of  the  Surface  due  to  Large 
Scale  Gravity  Waves 
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The  results  Indicated  in  Figures  3.2-7  to  3.2-11  are  encouraging 
and  indicate  that  the  Bragg  scattering  model  superimposed  on  a tilting  surface 
as  described  above  compare  favorably  to  the  TAGSEA  measurements.  Further, 

Pierson  and  Stacy's  model  was  derived  by  a compltely  independent  procedure, 
the  favorable  comparison  to  the  TAGSEA  results  tends  to  validate  both  sources. 

It  should  be  noted  that  Pierson  and  Stacy  compared  their  model  to 
NRL  backscatter  measurements.  They  were  able  to  get  very  good  agreement  except 
for  a factor  of  about  12  or  4ti  bias  between  their  predicted  backscatter  based 
on  their  model,  and  NRL  measurements.  They  implied  that  the  NRL  measurements 
were  too  low  by  a factor  of  12.  Since  NRL  data  represents  the  most  comprehensive 
set  of  clutter  measurements,  all  TSC  clutter  models  published  previously  have 
been  heavily  biased  by  the  NRL  data. 

3.2.2  Implications  of  the  Theoretical  Model  on  the  Statistical 

Fluctuation  of  Sea  Backscatter 

The  model  discussed  above  leads  fairly  naturally  to  a space  and  time 
varying  Rayleigh  fluctuation  model  as  proposed  by  several  authors  (e.g., 

Guinard  and  Daley  [20],  Sodergren  [8],  Tong  [2l],  and  Trunk  [22]).  The 
instantaneous  measured  cross  section  consists  of  a short  term  fluctuation 
which  is  determined  by  the  scattering  from  the  capillary  waves.  The  mean  level 
of  this  short  term  scattering  is  given  by  eq.  (6).  Note  that  the  expression 
depends  on  the  local  angle  of  incidence  and  the  mean  wind  velocity.  As  different 
parts  of  sea  are  examined,  the  local  angle  of  incidence  will  vary  as  the  gravity 
waves  and  swell  cause  tilting  of  the  surface.  Further,  spatial  variations  in 
wind  velocity  (turbulence)  cause  the  height  of  the  capillary  waves  to  vary 
accordingly.  As  noted  by  Pierson  and  Stacy  [5],  the  capillary  waves  build  up 
very  quickly  and  act  as  almost  instantaneous  tracers  of  the  local  wind  speed. 

Qualitatively,  such  a model  compares  favorably  with  the  TAGSEA  model 
proposal  in  Cl].  In  that  reference,  the  fast  fluctuation  component  was 
assumed  independent  from  FFT  frame  to  FFT  frame.  Thus  in  that  model,  the 
fast  component  decorrelates  in  9 milliseconds  or  less.  The  slow  fluctuation 
component  was  modeled  as  having  a spatial  correlation  of  300  to  1000  feet  and  a 
temporal  correlation  of  3 to  10  seconds. 
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Quantitative  comparisons  between  a theoretical  model  based  on 
Bragg  ucattering  and  Pierson  and  Stacy's  sea  surface  model  are  difficult 
to  make.  Generally  very  careful  interpretation  of  the  sea  surface  model 
and  Its  interaction  with  the  radar  system  is  required  to  obtain  a good 
quantitive  estimate.  Further,  formulas  resulting  from  the  model  usually 
require  fairly  involved  numerical  evaluations.  For  these  reasons,  no  attempt 
to  make  detailed  computations  will  be  made;  however,  some  relatively  quick 
calculations  can  be  considered, 


An  estimate  of  the  rate  of  fluctuation  for  the  fast  component  can 
be  obtained  by  considering  the  "phased  array"  model  for  the  diffuse  component  as 
discussed  earlier.  At  any  inBtant  of  time  the  diffuse  spatial  scattering 
appears  as  random  sidelobes,  produced  by  the  random  phase  errors,  caused  by 
the  capillary  waves  on  the  surface  as  depicted  in  Fig.  3.2-13. 


The  scattering  pattern  will  change  as  the  structure  of  the  capillary 
waves  changes,  in  addition  as  the  aircraft  flieB  along  it  samples  the  scattering 
pattern  in  space  at  an  angular  rate: 


• v 

n - - (8) 

For  a given  resolution  cell  size,  the  average  angular  width  of  a lobe  in  the 
scattering  pattern  can  be  expected  to  be 

n ■ | (9) 

TtiuB  if  the  surface  is  "frozen"  and  the  main  effect  is  the  aircraft  motion 
through  the  scattering  pattern,  the  correlation  time  is  approximately: 


T 


n 

n 


XR 

DV 


For  D»30m,  R“150Qm,  \*.03m,  V-130m,  the  correlation  time  is: 


T ■ .011  sec. 


(10) 
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For  this  case,  the  fast  fluctuation  component  decorrelates  sample* 
to-sample  in  agreement  with  the  TAGSEA  model.  Note  that  this  relation  will 
always  hold  if  the  azimuth  resolution  is  obtained  by  synthetic  aperture 
techniques.  However,  if  the  A/C  were  flying  towards  the  resolution  cell, 

the  significant,  angular  rate  would  be  in  the  vertical  plane,  and  would  be 


A V sin  cp 
* 


(ii) 


where  cp  is  the  grazing  angle.  Continuing  the  example, the  aircraft's  contribution  to 

correlation  time  is  ^ 

..."  T ■ rr-n — 7 — 7 ■ .06  sec.  maximum* 

D V sin  V 


,r*  Q 

at  10  grazing  angle.  Thus  for  the  forward  looking  cases,  the  correlation 
time  of  the  fast  fluctuating  component  might  be  significantly  longe-  than  for 
the  sidelooking  case  if  it  .were  not  for  wave  motion.  Note  that  the  theoretical 
model  would  indicate  that  in  nonsynthetic  aperture  cases  and  for  high  processing 
frame  rates,  the  fast  fluctuation  component  would  be  correlated.  The  TAGSEA 
model  Implies  (but  does  not  state  specifically)  that  the  fast  fluctuating 
component  is  independent  under  ail  conditions. 


The  theoretical  model  can  also  be  uaed  to  approximate  the  characteristics 
of  the  slowly  fluctuating  component.  Recall  from  eq.  (6)  that  the  mean  cross 
section  is  a function  of  local  wind  speed  and  local  grazing  angle.  Using 
Pierson  and  Stacy' a [5]  model  for  gravity  waves  it  is  possible  to  approximate 
the  variation  in  local  grazing  angle.  Further,  using  models  of  atmospheric 
turbulence  (e.g.,  Kalbaugh  [23]),  it  1b  possible  to  approximate  the  variations 
expected  due  to  these  variations  in  local  wind  conditions.  As  noted  earlier, 
detailed  computations  will  not  be  attempted  here.  However,  it  is  interesting 
to  note  that  if  the  rms  variation  in  the  cross  section  is  computed  using 
Eq.  (6)  and  averaged  over  a Gaussian  slope  distribution,  the  values  are  close 
to  0.2  as  shown  in  Table  111  for  a typical  case.  This  value  corresponds  to 
the  extreme  mean  variation  case  noted  in  the  TAGSEA  model. 


A final  comparison  between  the  TAGSEA  model  and  the  theoretical 
model  ia  in  the  underlying  distribution  for  the  rapidly  varying  component. 
The  TAGSEA  model  proposes  a slightly  nonrayleigh  distribution  to  model  the 
rapidly  fluctuating  component.  From  the  previous  discussion,  it  might  be 

* NOTE:  This  ignores  the  wave  dynamics  of  the  sea. 
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Table  III.  RMS  Variation  In  the  Long  Term 
Fluctuation 

(Theoretical  Computation) 


o Wind  Velocity  15  kt 
o Vertical  Polarization 


Grazing  Angle 


10 

20 

30 

40 

50 


RMS  Variation 
(dimensionless) 

.57 

.28 

.21 

.23 

.33 
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inferred  Chat  Che  theoretical  model  would  imply  a Rayleigh  distribution  for 
Che  rapid  fluctuations.  It  should  be  noted  that  the  Bragg  scattering  model 
noC  8 many,  equal,  independent  scatterer  model.  The  received  field  la 
modeled  as  the  spatial  Fourier  transform  of  the  capillary  wave  heights  over 
the  Illuminated  range- azimuth  resolution  cell.  In  a mathematical  sense,  the 
Fourier  transform  will  be  Gaussian  (Rayleigh  amplitude)  only  when  the  capillary 
waves  can  be  represented  as  a stationary  Gauasian  process.  As  a practical 
matter,  the  exact  distribution  of  the  capillary  waves  is  not  as  important 
as  the  assumption  of  statlonarity,  As  noted  above,  gravity  waves  and  swell 
can  produce  significant  nonstetlonarlty  by  tilting  the  local  eurfeca.  It  ie 
reasonable  to  ask  under  what  conditions  will  the  sampled  0°  fluctuations 
approximate  t Rsyloigh  (power)  distribution. 


Two  extreme  cases  are  of  interest.  In  the  first  extrema,  the 
resolution  cell  is  so  broad  that  the  tilt  and  gusts  are  averaged  out  and 
the  residual  oq  variations  are  negligible.  This  is  the  ueusl  low  resolution 
radar  limit.  At  Cha  other  extreme,  the  resolution  cell  ie  kept  fixed  relative 
to  the  surface,  and  tha  samples  are  taken  in  so  short  e time  that  the  sea 
surfaca  and  wind  gusts  ara  "frozen".  By  sampling  the  resolution  cell  from 
different  aspects,  one  can  obtain  independent  samples  from  tha  Rayleigh  dis- 
tribution which  correspond  to  that  cell.  Note  however,  that  tha  maximum 
number  of  samples  which  can  be  obtained  la  limited.  In  the  TAGSEA  case,  it 
was  shown  earlier  that  tha  maximum  independent  sampling  rate  is  about  .011  sec. 
For  e resolution  cell  of  size  30  m,  the  dominant  tilting  effect  is  from  gravity 
waves  of  approximately  the  same  size.  Thus,  the  tilting  should  occur  with  a 
nominal  frequency  of 


JO.  m M 

2tt  2Tt 


.22  He. 


The  available  number  of  samples  is  on  the  order  of  400,  end  it  is  difficult 
to  obtain  accurate  statistics  from  such  a limited  sample  size.  One  possible 
approach  would  be  to  normalise  the  statistics  to  the  mean  leval  of  the  400 
samples,  and  then  cumulate  statistics  on  the  normalized  date,  The  result  should 
approach  a Rayleigh  distribution.  It  should  be  noted  that  the  above  procedure 
does  not  correspond  to  the  N-type  normalization  used  in  reference  [l]. 
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Although  tho  data  was  normalized  by  range  gate  and  over  short  time 
periods,  it  should  ba  noted  that  for  histogram  analysis,  all  Doppler  bins 
were  cumulated  without  normalization.  The  total  Doppler  spread  corresponds 
to  approximately  3200  feet  along  the  aircraft  track;  thus  the  data  is  taken 
such  that  several  different  gravity  wave  systems  are  lumped  together.  The 
result  is  that  tho  N-type  fast  fluctuations  soon  by  the  TAGSEA  data  analysis 
correspond  to  a non-Rayleigh  distribution. 

It  is  quite  likely  that  the  differences  mentioned  above  are  largely 
academic  since  most  CFAR  systems  average  over  several  range  cells,  and  often 
over  several  azimuth  cells.  Thus  the  distribution  seen  by  the  CFAR  circuit 
will  be  non-Rayleigh  as  proposed  tv-  the  TAGSEA  model. 

In  summary,  the  TAGSEA  model  agrees  well  qualitatively  with  a time 
end  space  varying  Rayleigh  distribution.  The  theoretical  model  is  based  on  a 
composite  sea  surface  model  for  which  the  dominate  scattering  mechanism  is 
Bragg  scattering  from  the  capillary  waves.  Gravity  waves  and  swell  serve  to 
Introduce  nonstationarity  into  the  data  by  tilting  the  local  surface  in  a 
resolution  cell.  Additionally,  preliminary  calculations  indicate  significant 
quantitative  agreement  exists  between  the  theoretical  model  and  the  TAGSEA 
data  for  vertically  polarized  cases.  Since  the  sea  surface  spectrum  proposed 
by  Pierson  and  Stacy  [5]  is  brsed  on  direct  water-wave  measurements  rather 
than  radar  measurements,  the  quantitative  agreement  tends  to  lend  credence 
to  both  results.  Quantitative  results  for  horizontal  polarization  are  not 
as  good  at  loast  for  reflectivity, 

3.2.3  Comparison  of  the  TAGSEA  Fluctuation  Model  with  Other  Published 

Results 

Although  there  are  several  sources  for  published  results  on  the 
average  sea  clutter  cross  section,  there  are  fewer  published  results  on  the 
fluctuation  statistics  of  Oq.  Most  investigators  consider  only  the  composite 
distribution  (Type-A  or  second  order  model  in  the  nomenclature  of  [l]).  In 
a following  subsection,  the  TAGSEA  model  for  "Type-A"  distributions  will  be 
compared  to  several  other  experimental  results. 
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To  this  writer's  knowledge,  four  authora  ([24],  [22],  [2l],  and 
[8])  giv«  detailed  charactarlatica  of  a modal  which  attempts  to  daacrlba  the 
underlying  long  tarn  fluctuation  charactarlatica  aa  waa  dona  in  the  TAGSEA 
third  order  modal. 

Wright  [24]  conaldara  only  the  naan  croaa  aectlon  aa  waa  dona  In 
the  previous  section.  He  obtained  comparable  results  In  that  his  results 
for  vertical  polarization  agree  well  with  backscattar  data.  Trunk  [22] 
computed  the  statistics  of  the  composite  distribution  using  nearly  the  same 
formulas  as  were  given  in  the  previous  section.  He  alao  averaged  reaults  over 
a one  dimensional  slope  distribution,  but  baaed  his  calculations  on  the 
Kitaigorodskii  model  of  the  water-wave  spatial  spectrum^  Trunk's  reaults  for 
vertical  polarization  were  compared  to  the  TAGSEA  model  in  reference  [l]. 

There  it  was  shown  that  Trunk's  reaults  compere  well  with  the  TAGSEA  date. 

Tong  [21]  indirectly  inferred  the  characteristics  of  the  underlying 
distribution  by  matching  theoretical  calculations  to  observed  statistical 
distribution  functions. 

Sodergren  [8]  presents  experimental  reaults  on  both  the  composite 
distribution  (That  is,  the  distribution  which  reflects  both  the  long  and 
short  term  trends.)  as  well  as  giving  characteristics  of  the  long  term 
fluctuations.  Since  theoretical  calculation*  war*  presented  previouely,  only 
Sodergren' a date  will  be  compared  to  the  flow  fluctuation  charactarlatica 
of  the  TAGSEA  third  order  model. 

3. 2. 3.1  Comparison  to  Experimantally  Measured  Composite  (Type-A) 

Distributions 

In  this  subaaction,  the  TAGSEA  model  la  comparad  to  the  fluctua- 
tion atatlatics  measured  by  other  experlmentors . The  date  to  be  comparad 
corresponds  to  the  TAGSEA  second  order  modal,  that  is,  the  compoeita  distribu- 
tion arising  from  both  the  long  term  end  the  short  term  fluctuations. 

Although  several  raferancea  art  availabla  for  obtaining  tha 
fluctuation  atatlatics,  moat  raferancas  are  not  directly  comparable  to  the 
TAGSEA  maaauramenta.  In  particular,  Brooks  and  Brooks  [25]  show  that  the 

Sea  Pierson  end  Stacy  [5]  or  Trunk  [22]  for  a description  of  tha 
Kitaigorodskii  water  wave  apactrum. 
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published  reference*  exhibit  e moderately  strong  trend  of  Increasing  skewness 

In  the  composite  distribution  versus  gracing  angle.  Figure  3.2-14  taken  from 

C29]  shows  that  below  about  8 or  9 degrees  grazing  angle,  the  ratio  of  the  84 

_k 

percentile  to  the  50  percentile  increases  as  sin  v|r . Thus  data  taken  at  low 
grazing  angles  should  exhibit  significantly  higher  tails  than  the  TAGSEA  data 
which  was  taken  at  grazing  angles  of  5°  and  greater.  This  la  demonstrated  in 
Tables  IV,  V,  and  VI.  In  those  tables,  the  TAGSEA  model  is  compared  with 
data  taken  by  Bishop  [10]  and  Rivers  [11] ■ Both  seta  of  data  were  taken  at  low 
grazing  angles  in  shallow  water.  On  Figure  3.2-14,  Bishop's  data  yields  high 
values  of  O'  (7-17  dB)  and  as  seen  from  Tables  V and  VI,  the  distributions 
lie  significantly  above  the  TAGSEA  model.  Rivers  data  on  the  other  hand 
generates  relatively  low  valuee  of  o'  (4-6  dB)  but  still  lies  significantly 
above  the  TAGSEA  model. 

The  above  comparisons,  obviously,  provide  no  information  on  the. 
validity  of  the  TAGSEA  model,  lhey  do,  however,  serve  to  emphasize  the 
limitations  of  the  model  as  stated  in  [l],  The  model  was  derived  to  apply 
only  to  the  conditions  for  which  the  graateat  amount  of  data  was  available. 

Thus  the  TAGSEA  model  applies  only  to  moderate  grazing  angles  (3‘-43*),  high 
sea  state  (4-3),  and  deep  water  data.  As  the  above  comparisons  indicate, 
conditions  which  deviate  significantly  from  the  TAGSEA  experiment  are  likely 
to  require  significantly  different  models.  As  reference  [23]  indicates,  the 
greatest  affect  may  be  grazing  angle. 

A review  of  the  available  references  indicates  that  references 
[6]-[9]  and  [26],  [27]  have  data  for  grazing  angles  above  5 degrees.  From 
this  data  base,  results  based  on  the  NRL  four  frequency  radar  must  be 
eliminated  since  that  radar  la  a low  resolution  radar.  Although  the  pulse 
length  la  reasonably  small  for  that  radar,  the  beamwldth  is  quite  wide  (*■  5°) 
and  the  azimuth  resolution  call  is  broad.  Upon  inspection  of  the  references, 
it  was  found  that  many  of  the  references  presented  the  same  date.  Thus  after 
eliminating  the  redundancy,  it  was  found  that  only  Trunk  and  George  [6], 

Schmidt  [7],  Sodargren  [a] , and  Boeing  [9]  represent  nonredundant,  high- 
resolution,  moderate-grazing  angle  data  for  comparison  with  the  TAGSEA  model. 
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Grazing  Angle  (deg.) 


Table  IV.  Compariaon  of  the  TAGSEA 
Modal  ulth  Bishop  [lO] 


o Vortical  Polarization 

o Level  in  dB  re  mean 

o Bishop  data  at  1°  grazing 
45  x 40  m 


Prob.  of 
Excaad. 

TAGSEA 

Typical 

TAGSEA 

Extrema 

Biahoo 

Data 

SS-4-5 

SS-2-3 

SS-3 

SS-6 

io"1 

3.6 

3.6 

4.0 

4.5 

4.0 

4.5 

10‘2 

7.1 

7.0 

9.0 

10.0 

9.0 

10.5 

io"3 

9.3 

9.4 

11.5 

13.0 

13.6 

15.5 

-4 

10  4 

10.9 

11.5 

13.8 

14.6 

16.7 

19.7 

io"5 

12.2 

13.4 

15.0 

16.0 

19.5 

21.3 

IO’6 

13.3 

15.2 
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Table  VI.  Comparison  of  the  TAGSEA  Model  with 
Rivers  C 11 3 


o Grazing  Angle  « 1.5° 
o Echo  Area  **  12  m x 30  m 
o Level  in  dB  re  mean 
o Wind  Speed  9.5  kts  (SS  3) 


0 

X-band, 

Upwind 

Vertical  Polarization 

Prob.  of 
Exceed. 

TAGSEA 

Typical 

TAGSEA 

Extreme 

Rivers 

TAGSEA 

Tvoical 

TAGSEA 

Extreme 

Rivers 

lO"1 

3.5 

3.4 

3.6 

3.6 

3.6 

3.8 

10-2 

7.6 

8.1 

9.8 

7.1 

7.0 

7.8 

10"3 

10.4 

11.6 

14.6 

9.3 

9.4 

11.1 

i 

j 

10-4 

12.6 

14.5 

19.1 

10.9 

11.5 

14.0 

K-* 

o 

1 

Oi 

14.6 

17.1 

12.2 

13.4 

10'6 

16.4 

18.5 

13.3 

15.2 

■ 

i 

[ 

f 
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The  chrracteristics  of  the  data  which  can  be  compared  to  the 
TAGSEA  data  are  given  in  Table  VII.  Of  these  references,  Sodergren's 
results  were  already  compared  with  the  TAGSEA  data  in  [l].  As  can  be 
seen  from  Figure  3.2-15  taken  from  [l],  his  reoults  compare  well  with 
the  TAGSEA  second  order  uiodel.  A more  complete  discussion  of  Sodergren's 
model  versus  the  TAGSEA  model  is  given  later. 

The  Boeing  results  are  preliminary,  have  only  been  reported  in- 
formally, and  are  at  S-band.  For  these  reasons,  it  would  be  premature  to 
report  on  them  in  detail.  However  they  did  observe  non-Rayleigh  distributions. 
When  they  fit  the  linear  envelope  with  a Weibull  distribution  they  obtained 
an,  "v"  parameter  of  *1.75.  This  result  is  not  inconsistent  with  the  TAGSEA 
results  which  obtained  a Weibull  parameter  of  *®.8  for  the  Bquare  law  de- 
tected output  [4J.*  A linear  law  Weibull  parameter  of  1.75  corresponds  to  a 
square  law  Weibull  parameter  of  1.75/2  - .875. 

Table  VIXI  compares  distributions  given  by  Trunk  and  George  [6] 
and  Schmidt  [7]  with  the  TAGSEA  model . As  can  be  seen,  the  data  in  Table  VIII 
appears  to  be  significantly  higher  than  the  TAGSEA  model.  Trunk  and  George's 
results  appear  to  be  quite  high,  ar.i  in  fact  may  have  been  selected  for  that 
reason  since  they  were  reported  in  an  Introduction  to  a paper  on  detection  in 
nongaussian  clutter.  Another  jossibillty  is  that  the  extromely  short  pulse 
length  (3  m)  of  that  data  may  '-untribute  to  the  high  tails.  Schmidt’s  data, 
taken  at  higher  sea  state  with  the  same  radar  does  not  have  as  high  tails  as 
the  data  presented  by  Trunk  and  George.  Note  that  Schmidt's  data  contradicts 


* The  approximate  1 ilue  for  N corresponding  to  the  Weibull  distribution  can 
be  obtained  from  the  statistical  parameters  given  in  Section  1.2  of  ref. 
[4j  . For  a given  value  of  Q,  the  percentile,  X,  corresponding  to  that 
value  is  given  by: 

• x _ U(-In-Q)l/N 

r<i  + n*1) 

where  N 1b  the  Weibull  parameter,  p is  the  mean,  and  T is  the  gamma 
function.  In  Section  1.2  of  [4],  the  data  is  normalized  so  that  p-10. 
TIiub  for  N«.8  and  Q-10"^,  then  X*22.7  dB.  This  corresponds  to  row  #14 
of  the  table  given  in  Section  1.2.1  of  [4],  in  that  reference  the  values 
vary  from  about  21.5  dB  to  28  dB, 
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Figure  3.2-15.  Comparison  of  the  TAGSEA  Composite 
Distribution  to  Sodergren's  Data  IB]. 
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Table  VIII.  Comparison  of  the  TAGSEA  Model  with  Trunk  and  George 
[6]  and  Schmidt  [7] 


Vertical  Polarization 

Horizontal  Polarization 

Prob. 

of 

Exceed. 

TAGSEA 

Typical 

TAGSEA 

Extreme 

Trunk  & George 
[6] 

SS-1  SS-2-3 

Schmidt 

[7] 

SS-5-6 

TAGSEA 

Typical 

TAGSEA 

Extreme 

Schmidt 

[7] 

SS-5-6 

•'r  ~ 
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Trunk  and  George's  data  which  shows  higher  tails  at  higher  Bea  states. 
Further,  the  trend  of  higher  tails  at  higher  sea  states  is  not  completely 
consistent  with  Bishop's  vertical  polarization  data  (Table  IV). 

However  Bishop  also  shows  atypical  cases  which  contradict  his  trend  of 
lower  tails  with  higher  sea  state. 


It  does  not  appear  possible  to  prove  or  disprove  the  validity  i 1 

of  the  TAGSEA  model  on  the  basis  of  other  experimental  data.  It  is  possible  ! 1 

to  find  data  which  agrees  and  data  which  disagrees  with  the  model.  When 
the  rather  large  list  of  references  which  give  distribution  functions  is'  , ■ 

sifted  to  obtain  those  data  which  can  be  directly  compared  to  the  TAGSEA  .■ 

model,  it  is  found  that  only  four  data  seta  are  left.  Of  these  data,  \ 

Sodergren's  data  and  the.  Boeing  data  have  the  strongest  similarity  to  the  ' 

TAGSEA  modei.  The  other  data,  taken  by  the  NRL  20  nsec  radar  have  higher  j 

tuili  for  the  composite  distribution.  Thus  it  would  appear  to  be  a conscr- 
vative  design  practice  to  at .least  consider  the  effect  on  system  performance 
of  bails  which  are  2-3  dB  higher  at  the  10  ^ point  on  the  distribution. 

3. 2.3. 2 Comparison  with  Sodergren's  Fluctuation  Model 

Ab  noted  earlier,  Sodergren  is  the  only  source  for  an  experimentally 
determined  model  comparable  to  the  TAGSEA  third  order  model.  As  has  been 
noted,  Sodergi-en's  distribution  functions  compare  favorably  with  the  TAGSEA 
second  order  model.  In  this  section,  comparisons  related  specifically  to 
the  third  order  models  will  be  made. 

Sodergren,  TAGSEA  and  several  theoretical  models  propose  a composite 
fluctuation  model  which  represents  the  cross  section  function  statistics  as 
the  product  of  a slow  and  a fast  fluctuation.  That  is, 

°o  * °0  * S ’ F (12)  | 

The  fast  fluctuation,  F,  represents  the  "instantaneous"  Bragg  scattering  from 

the  capillary  waves,  while  the  slow  fluctuation,  F,  represents  the  variations  in  j 

the  average  backacattar  strength  of  the  capillary  waves  due  to  local  surface 

tilting  and  wind  gusts.  r 

i 
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Table  IX.  Compar’aon  of  the  TAGSEA  Third 
Order  Model  with  Sodergren'a  Model 


Rapid  Fluctuation 

0 

Statistics 

0 

Correlation 

Slow  Fluctuation: 

0 

Statistics 

0 

Correlation 

o Space 

o Time 

TAGSEA 

Sodergren 

Slightly  Non-Rayleigh 

<j 

Rayleigh 

Independent 

Dependent  on  SS  and  X 

Normal 

i 

1 

1 

Log  Uniform  or  Log  ! 

Normal  ,j 

Exponential 

\ 

Constant  • j 

Exponential 

Constant  j 

Table  IX  gives  a qualitative  comparison  between  the  TAGSEA  modal 
and  Sodergren's  model.  As  indicated  by  the  table,  the  details  of  the  two 
models  differ  in  that  the  assumed  distribution  functions  and  correlation 
functions  are  different.  However,  quantitative  comparisons  should  be  made 
before  any  conclusions  can  be  drawn, 

First  consider  the  slow  fluctuation  model  correlation  functions. 
Neither  model  is  definitive  in  this  respect.  Sodergren's  model  assumes  short 
processing  segments  and  times  which  are  localised  in  space;  thus  the  long  term, 
large  Bpatial,  variations  of  the  fluctuations  were  not  considered.  TAGSEA 
gives  a tentative  model  based  on  a relatively  small  amount  of  supporting  data. 
The  model  assumes  temporal  correlation  times  of  3-10  seconds  and  spatial 
correlation  distances  of  90  to  300  tn.  High  resolution  radar  sets  with 
resolution  cells  on  the  order  of  IS  to  30  m which  perform  requitted  CFAR 
averaging  on  10  to  20  cells  should  not  be  greatly  sensitive  to  either  modal. 
Radar  sets  which  average  over  larger  areas  should  be  evaluated  against  the 
TAGSEA  model  since  it  provides  a more  difficult  test  thavi  the  Sodergren  model. 

Next,  consider  the  statistics  of  the  long  term  fluctuations.  For 
Bmall  variance  of  the  long  term  fluctuation,  the  normal,  log  normal,  and  log 
uniform  distributions  do  not  differ  significantly.  Sodergren  shows  that  tha 
relation  between  the  standard  deviation  of  the  slow  fluctuations,  a^,  and 
the  composite  standard-deviation-to-mean  ratio,  SDMR,  is  given  by 

O - Cl/2  (SDMR2  - l)]l/,Z 

8 

for  a time  varying  Rayleigh  distribution.  This  result  does  not  depend  on  the 
distribution  of  S.  Thus  for  Os  ■ .1  and  .2  one  gets 

SDMR  - 1.01  TAGSEA  typical 

- 1.04  TAGSEA  extreme 

Sodergren  gives  SDMR  versus  grazing  angle  and  sea  state.  This  la 
given  as  Figure  3.2-16  and  the  TAGSEA  SDMR  as  computed  above  are  marked.  This 
figure  indicates  that  the  TAGSEA  model  would  produce  a low  value  for  the 
fluctuation  statistics.  Note,  however,  such  a comparison  is  not  accurate. 
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Figure  3 . 2-16 . Adjusted  Standard  Deviation- to-Mean  Ratio 
(From  Sodergren  [8]) 


As  an  example,  for  run  702,  which  corresponds  Co  Che  extreme  TAGSEA  model 
for  horizontal  polarization,  the  SDMR  computed  directly  from  the  data  la 

SDMRlQ2  - 1.405 

The  apparent  contradiction  lies  in  the  fact  that  the  TAGSEA  model 
is  a time  varying  non-Rayleigh  diotribution.  Thus  the  TAGSEA  model  takes  most 
of  the  fluctuation  out  of  the  slow  fluctuation  portion  and  puts  it  into  the 
rapid  fluctuation.  As  noted  earlier  in  the  theoretical  diacusaion,  this  is 
primarily  a philosophical  point.  In  any  practical  system,  the  extreme  difficulty 
of  drawing  samples  from  the  samo  resolution  cell  under  identical  conditions 
makes  it  nearly  impossible  to  physically  observe  the  underlying  Rayleigh 
distribution.  Thus  correctly  using  either  model  should  produce  equivalent  CFAR 
designs  and  performance  valuations. 

Note  that  tha  correct  application  of  a time  varying  Rayleigh  model 
requires  sampling  from  slightly  different  resolution  cells  as  usually  happens 
in  actual  data  acquisition.  This  is  a subtle  point  in  the  application  of  the 
time  varying  Rayleigh  distribution  which  can  be  overlooked.  The  TAGSEA  model 
is  somewhat  preferred  since  it  has  already  included  some  of  this  fluctuation 
into  the  rapid  fluctuation  model. 

The  rate  of  rapid  fluctuations  is  different  between  Sodergren's 
model  and  the  TAGSEA  model.  The  TAGSEA  model  assume  "white"  rapid  fluctuations 
while  Sodergren's  model  assumes  the  bandwidth  of  the  fluctuations  are  given  by 

BW  - ^ 

where  \ is  rf  wavelength  and  SS  is  hydrographic  sea  state.  M sea  state  5, 
and  X band,  Sodergren's  model  gives 

BW  - 250  Hz 

At  the  TAGSEA  FFT  frame  rate  (*>  100  Hz),  this  corresponds  to  independent 
samples  frame-to-framo.  Thus  within  the  measurement  limitations  both  models 
agree. 
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1C  should  b«  notsd  that  from  eha  thaoratlcal  discussion  given 
earlier,  for  high  speed  radar  platforms,  the  bandwidth  of  the  fluctuations 
should  be  govarnod  by  wave  motion,  the  radar  resolution  cell  and  the  motion 
of  the  platform  relative  to  that  cell.  This  is  in  contradiction  to  Sodergren's 
model  (which  assumes  that  the  bandwidth  is  governed  by  sea  state)  and  Che  TAGSEA 
modal  which  assumes  che  rapid  fluctuations  are  white.  The  latter  assumption 
is  safe  for  low-data-rata  systems,  but  for  high-data-rete  systems  the  shape  of 
the  fast  component,  autocorrelation  function  (or  spectrum)  may  be  important. 

For  both  situations,  the  slow  component  can  contribute  appreciable  to  the 
distribution  statistics. 


As  a final  point,  it  should  be  noted  that  the  TAGSEA  fluctuation 
s»del  is  lndepandent  of  sai  state  and  gracing  angle.  This  la  conslstant  with  the 
TAGSEA  data  as  illustratad  in  'iablj  X.  Thare  tha  data  sata  are  grouped  such 
that  only  saa  stata  varies  within  e group.  Note  thet  generally  (except  for 
runs  401  and  403),  the  effect  of  sea  state  and  gracing  angle  is  small.  Thus 
Chs  TAGSEA  modal  in  conslstant  relative  to  its  data  sat.  Sodargrsn  shows 
that  the  standard-devlatlon-to-maan  ratio  varies  ae  a function  of  gracing  angle 
end  eee  atate  (see  Fig.  3.2-16).  However,  for  gracing  angles  above  20°,  the 
effect  is  not  significant.  At  low  gracing  angles  both  Sodergren's  date  and 
Biahop's  data  exhibit  fairly  strong  saa  state  affects.  Thus  the  TAGSEA  model 
should  only  be  applied  at  moderate  grazing  angle  and  high  saa  state. 
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